double precision function terfc(arg)
    implicit none
    double precision arg
    terfc = 1 - erf(arg)
end function terfc
